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1. INTRODUCTION 

Smart antennas (SA) were called also adaptive arrays [1, 2]. This term represents the capacity of the 
antenna to adjust to the condition of the communication channel in which it works. SA arrays with adaptive 
beam-forming capability are very efficient in suppressing interference and multi-path signals and improving 
signal quality by reducing the fading effects [3]. SA focuses on the notion of spatial processing, which is the 
key concept in adaptive antenna systems [4]. Arrays may produce the desired radiation characteristics by 
properly exciting each element with certain amplitudes and phases to maximize the signal from the desired 
users [5]. This paper analyzed various array geometries such as uniform linear (LA), uniform rectangular 
(PA), and circular arrays (CA). The LA has excellent direct and narrowest main beam lobe in a wanted 
direction, but in all azimuthal directions it does not work equally well, a major drawback of the PA is the 
question of presence on the opposite side of an additional large lobe of the same strength [6] the symmetry of 
the CA structure provides an obvious benefit since, it has no edge components, without a major change in the 
beamform, directional patterns synthesized with a CA can be rotated electronically in the array plane [7]. 
Many algorithms have been introduced to applicants on SA [8, 9], one of the most public algorithms is least 
mean square (LMS) A commonly employed algorithm due to its low computational complexity and ease of 
implementation [10]. Other less square algorithms, such as recursive least squares (RLS), conjugate gradient 
(CG), will converge more easily and have a lower stable mean square error (MSE) than LMS [11]. However, 
their high computational complexity makes them inadequate for many real-time applications [12, 13]. For 
adaptive filtering applications, the RDES algorithm is investigated called the Euclidean direction search 
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(EDS) since it searches cyclically for the minimum along one Euclidean direction at a time [14, 15]. These 
algorithms are investigated for decision feedback equalization in a mobile communication channel with 
additive noise and multipath. The RDES algorithm has a slower convergence rate than the Euclidean 
direction search (EDS) but a lot faster than the LMS [16, 17]. Also, the performance in terms of bit error rate 
(BER) and convergence rate is comparable to those of the RLS algorithm but at a much lower computational 
cost. The achievement of this work is organized by using the RDES algorithm [18, 19], with different array 
geometry for mobile communications. The result is analyzing by the RDES algorithm implemented by an 
EDS procedure to different array geometry and study the mathematical effect of array radiation pattern with 
applying AWGN and Rayleigh fading channels on interference suppression capability (null depth), mean 
square error (MSE), stability and convergence rate. 


2. RADIATION PATTERN FOR ANTENNA ARRAY SYSTEM 

The radiation pattern in most cases defined in the far field area and is interpreted as the directional 
coordinate's function [20]. Antennas with a given radiation pattern can be configured to create a various 
radiation pattern in a system with different dimensions. For an antenna array of identical elements, the 
radiation pattern of the array can be defined according to the pattern multiplication theorem [21]: 


Array pattern = Array element pattern x Array factor 


In which the pattern of an array element is the pattern of the single array element, and the array 
factor is a feature that depends only on the array geometry and the excitation of the elements (amplitude and 
phase). 


2.1. Linear array 

Figure 1(a) illustrated LA system, for simplification purposes assumption an equally spaced 
elements and equal amplitudes. However, each successor element has a (5) progressive step leading current 
excitation compared to the past one (6 is the phase in which the current leads the current of the preceding 
element in each element). Assuming that the nth element leads the (n-1) element by an electrical phase shift 
of 6 radians and can implement by shifting a phase of the antennas current for each individual element axis 
and the phase shift 6 is equal to zero. 

The AF for N elements can be considered as: 


AF = [1 + ej(kdsin 0+6) + ej2(kdsin 0+6) .. @ N-19 (kdsin 0+8)J]T (1) 


and, k = 2 * 2i is the phase constant, d is element spacing, O is arrival angle. In (2) can be expressed by, 


N 
AF = > Wee ov (2) 


n=1 
where, Y = kd sinO + 6 and 6 = —(kd + sin@) which is the progressive phase, k is the wavenumber and, 
Wn = a, * e? (3) 


2.2. Circular array 
Figure 1(b) shows CA in the x-y plane, with uniform angular distribution between elements of 


value Ø, = ae. The nth array element is positioned with the phase angle @, element at the radius (a). 


Every element can get a related weight w, and phase 6,. Also assuming far field conditions and conclude 
that the observation point such that the position vectors rand 7, are parallel. The AF can be found in a similar 
procedure as was calculated with the LA as [22] 


N N 
F= > We i kap P +n) a > wpe (kasin8 cos (p-pn)+ôn) (4) 


n=1 n=1 
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2.3. Planar array 
This type is used to perform 2D beamforming (in an azimuth and elevation angles) with horizontal 


spacing of d, and vertical spacing of d, is shown in Figure 1(c). The AF for PA can be expressed as 
combining the AF of two LAs [23]. Pattern multiplication can be used to calculate the pattern of the entire 
MXxN antenna array. Utilize pattern multiplication would have: 


M N 
AF = AF,. AF, = ` amel cose the) > peter D (kdysinðsing+£y) (5) 
m= n=1 
M N 
= 2 > We pete kdysinOcosp+px)+(n-1)(kdysinbsing+8y)] (6) 











(a) (b) (c) 


Figure 1. (a) Linear array, (b) Circular array, (c) Rectangular array 


3. ADAPTIVE ANTENNA ARRAY USING RDES ALGORITHM 

Figure 2 illustrated the block diagram of the array system. The signal S(k) and interferers 
1, (k), I2(k), ... In (k) are received by an array of N elements with N potential weights. Each received signal at 
elements n includes the noise. Time is represented by ky, time samples. The weighted array output can be 


given as: 


Desired direction 
(main beam) sik) 








Antenna 1 


Antenna 2 


Beamforming 









signal 
Antenna 4 y 
a 
direction (null) 
I{k} atk) 
EE 
algorithm ial = 
Figure 2. Block diagram for smart antenna 
y = Wik). x(k) (7) 
and, 
l(k) 
X(k) = āS(k) + [a, ā; .... ay]. I, o + (k) = Xs(k) + x,(k) + n(k) = input signal (8) 
I(k) 
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Xs(k) is the desired signal vector, X; (k) is the interfering signals vector. n(k) the zero-mean Gaussian noise 
for the channels. a,Are the N element array steering vector for O; is the arrival direction. The error signal 
is e(k) which defined as the difference of desired signal d(k) and output signal y(k) in the method of 
steepest-descent, on the other hand, a weight vector is chosen to minimize the ensemble average of the error 
squares [17]. The step-size parameter u governs the convergence of algorithm because it is directly 
proportional to it, therefore, imperative to choose a step-size in a range that ensures convergence [1, 24, 25]. 
The reformulation written with RDES, the error signal, and the weight vector can give as: 


e(k) = d(k) — W#(k)X(k) (9) 
wk + 1) = W(k) + pe*(k)X(k) (10) 


where the input and desired signals are X(k) = [x, Xz ....Xy]'&d = [d(1)d(2) ....d(K)]?. The error signal 
e(k) can be expressed as: 


M 
e(k) = d(k) - È wx a1) 


i=1 


Considering the samples k-L, k-L+1, k-L+2 ......k, where L>M, in (11) would write as: 


Elk) = d(k) — U(K)w(k) (12) 
where, 

d(k) = [d(k), d(k — 1),d(k — 2)... d(k -L+1)]" (13) 

U(k) = [ū; (k), 0) (k) ... dy (K)] (14) 

Tjk) = [x\(k), x(k — 1)... xj(k — L + J" (15) 


The new approximation error at time k is given by, 

E(k) = d(k) — U(k— 1)W(k - 1) (16) 
By update only one weight in W(k — 1) the new error can be written as [26]: 

E(k) = d(k) — [UKW — 1) + UK) WINS (k) Figo] (17) 


jo (k) is the index of the weight to be updated in the zeroth P-iteration iteration at time k and Fj, is Mx 1 
vector with 1 in position j and 0 in all other positions. After selecting index jo(k), the weight is given by, 


Tw k) = Wag (k — 1) + WN) (18) 


where Wj nS "(ka projection value of the error €)(k) and can be givens as: 


< &o(k)ū a(k) > 


gordate (k) = 
(k) = 2 
i līs] 


(19) 


when P-iteration is equal to zero the updates the array weight written with the step-size of as follows: 
= = — upd 
Wek) = Wk- 1) + proes Wty CX) Fpa (20) 


The optimized value of [gprs can calculate by applying the same parameter (p) used in (10). 
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4. RESULT AND DISCUSSION 
The simulation results are contrasted with three classical geometries. The number of elements is 
illustrated in Table 1 for comparison purposes in all simulations presented in this paper. 


Table 1. Initialization parameters 
No.elements Displacement distance 








Linear array 10 0.5A(x-axis) 
Circular array 10 2m/N 
Planar array 5*5 0.54 * 0.5A(x-y plane) 





An array consisting of isotropic elements and will be work under the following conditions: Input 


signal S(k) = sin(2tft(k)) with f = 1 = 900 MHz and t = (1: K) * T where K is represent the number of 
8 T K p 


the sample intervals and T is the period and that will be the desired signal d(k), (d(k) = S(k)), desired DOA 
is 6) = 0° and an interfering signal, 11 with DOA (0, = 45°). 

The initial step size parameter pọ = 1. Zero mean Gaussian noise with variance o2 = 0.001 would 
be added to the input signal for every element in the array system. The signal to noise ratio (SNR) and signal 
to interference ratio (SIR) are set at 30 dB and 10 dB respectively. The step size parameter is calculated 
according to (10) as: 


u = 1/(4 * real(trace(R,x))), where Ry, = X * X' 


an analytical model operating in the presence of two types of noise is discussed here to clarify the 
performance of the RDES algorithm under the effect of additive white Gaussian noise (AWGN) and 
Rayleigh fading channels. 


4.1. Additive white gaussian noise channel results 

The fundamental noise model used in information theory to imite the effect of several random 
processes that occur in nature. As apparent in Figures 3, Figures 4, and Figures 5 that LA and PA have more 
side lobe than CA but at the same time the null depth of LA and PA is larger than CA, also the sidelobe 
amplitude PA is less than LA and CA array configurations. 
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Figure 3. Linear and polar radiation patterns for linear array (LA), all geometry worked out with 50 
samples/AWGN 
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Figure 4. Linear and polar radiation patterns for circular array (CA), all geometry worked out with 50 
samples/AWGN 
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Figure 5. Linear and polar radiation patterns for planar array (PA), all geometry worked out with 50 
samples/AWGN 


The RDES has a rate of convergence faster by comparison with the previous algorithm like LMS which 
needs 50 iterations to start adaptation processes, NLMS algorithms that need more than 12 iterations, and 50 
iterations for CMA and so on. Figures 6(a), Figure 7(a), and 8(a) could illustrate the array output of PA and 
CA is a good estimate of the desired signal. Figures 6(b), Figure 7(b), and Figure 8(b) shows the magnitudes 
of the estimation complex weights of the, LA, CA and PA configurations, which is the magnitude of the 
weights for each element in array. 

The curves of the MSE are presents in Figures 6(c), Figure 7(c) Figure 8(c) that shows the 
convergence speed using PA geometry is faster than LA and CA configuration. Moreover, the array output 
using PA was a good estimate of the input signal compared with LA and CA but, at the same time, it can see 
that the error magnitude of CA is less than both LA and PA. The results show the performance of the RDES 
have so good stability and could notice the adaptation processes starts approximately after 5 samples only. 
The numerical result can be shown in Table 2. 
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Figure 6. Linear array related curves, (a) Signal tracking of desired and output curves, (b) Weight 
magnitudes, (c) MSE for linear array (LA) worked out with 50 sample/AWGN 
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Figure 7. Circular array related curves, (a) Signal tracking of desired and output curves (b) Weight 
magnitudes 
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Figure 7. Circular array related curves, (c) MSE for circular array (CA) worked out with 50 sample/AWGN 
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Figure 8.Plannar array 
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Table 2. Additive white Gaussian noise results 
Depth null (dB) SLL (dB) Weight stability Maximum mean square error 








LA -61 -13 8 iterations 1.4089e-04 
CA -60 -10 5 iterations 0.0055 
PA -96 -23.6 7 iterations 1.2625e-04 





4.2. Rayleigh fading channel results 

Rayleigh fading is the statistical model that effect on the radio signal of a propagation environment, 
is most applicable when several objects in an environment scatter the radio signal and that lead to no 
dominant propagation along a line of sight between the transmitter and receiver, such as that used by wireless 
devices [27]. The magnitude of a signal that has passed through such a transmission medium (a 
communication channel) will vary randomly is assumed. In the performance evaluation the following 
conditions are considered: The signals arriving at each antenna element, for the desired and interfering signal, 
subject to independent Rayleigh fading. Both the wanted signal and Interfering signal have the same 
amplitude as are emanating from 0° and 45°. Each simulation involves a run (data rate) of 500 Kbits. 
Doppler frequency fg of 117 Hz, corresponding to the mobility of 140 km/h at 900 MHz will be used in this 
simulation to present a worst case process. All antenna geometry is implemented when the number of 
samples is set at K=60. 

Figures 9, Figure 10, and Figure 11 show that the performance of the RDES which works better 
when compared with planar array than the LA&CA array. 
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Figure 9. Linear and polar radiation patterns for linear array (LA), all geometry worked out with 
60 sample/Rayleigh channel 
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Figure 10. Linear and polar radiation patterns for circular array (CA), all geometry worked out with 
60 sample/Rayleigh channel 
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Figure 11. Linear and polar radiation patterns for planar array (PA), all geometry worked out with 
60 sample/Rayleigh channel 


And as shown in Figures 12(a), Figures 13(a) and Figures 14(a) the convergence speed using LA is 
slower than CA configuration. Moreover, the array output using the CA algorithm is a good result of the 
input signal compared with the LA configuration. These figures confirm that the PA configuration has faster 
convergence speed and good output array estimation. . Figures 12(b), Figures 13(b), and Figures 14(b) shows 
the magnitudes of the complex weights of LA, CA, and PA configurations. From Figures 12(c), Figure 13(c), 
and Figure 14(c), the LA antenna converges after more than 30 iterations, CA antenna converges after 20 
iterations, whereas the PL antenna convergence can occur at about less than 10 iterations. 
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Figure 12. Linear array related curves, (a) Signal tracking of desired and output curves, 
(b) Weight magnitude 
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Figure 12. Linear array related curves, (c) MSE for linear array (LA), all geometry worked out with L=6, 
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Figure 13. Circular array related curves, (a) Signal tracking of desired and output curves, (b) Weight 
magnitudes, (c) MSE for circular array (CA), all geometry worked out with L=6, K=60 sample/Rayleigh 
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Figure 14. Plannar array related curves, (a) Signal tracking of desired and output curves, (b) Weight 
magnitudes, (c) MSE for plannar array (PA), all geometry worked out with L=6, K=60 sample/Rayleigh 
channel 


The results could present that the performance of the PA is better and faster than the other shapes. 
The numerical results with Rayleigh channel are shown in Table 3. 


Table 3. Rayleigh fading results 
Depth null(dB) SLL(dB) Weight stability Maximum mean square error 








LA -62 -13 40 2.4865e-04 
CA -48 -6 20 0.0010 
PA -90 -24 10 5.7605e-04 





Figure 9, Figure 10 and Figure 11 illustrates that the PA generates the deeper null which about 
-90 dB towards the interferer, the LA and CA configurations provide the deep null of about -62 dB and 
-48 dB at K=60. It is apparent from this array pattern in Figure 9, Figure 10, and Figure 11 the CA array 
provides the lowest array side lobes number while PA provides the lowest array SLL compared with the 
other configurations. Multipath and their effects on antenna array pattern with using RDES algorithm also 
have been introduced, Figure 15 illustrate linear and polar plot for antenna array pattern curve under 
multipath signals. 
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Figure 15. Linear & polar plots, for radiation patterns with 8 elements linear array where a-desire 
(S)=0/Interference (I)=45, -20, b-S=0, 35/I=60, -30, c-S=0, 30/1=50, -30, -70, d-S=0, -20, 30/I=50, -45, -70, 
worked out with 50 sample/AWGN 


As illustrated in Table 4, the effectiveness of the RDES algorithm appears to guidance or steer the 
mean pattern beam in the direction of users and steering the depth extinct to interference directions. 


Table 4. Analyze multipath effect 
Desired DOA(S) Interference DOA(I) Depth null (dB) 








a 0 45 &—20° -48.5&-59.2 

b  0°&35° 60 & — 30° -61&-41.7 

c 0°&30° 50°& — 30°& — 70° -38.1&-42.3&-62.7 
d 0 &30°&— 20° 50&—45 &—70 -36.7&-51.2&-55.1 





5. CONCLUSION 

The renewal in this paper is the use of the RDES algorithm in different antenna array geometries, 
also consideration the RDES algorithm's ability to operate under the effect of AWGN channel and Rayleigh 
fading environment. The result appears robust and remains stable when working in the case of Rayleigh 
fading and even when additive white Gaussian noise (AWGN) corrupts their reference signals. Via 
simulation of AWGN channel smart antennas the LA, PA, and CA begin to converge with iteration number 
of fewer than 10 samples. These configurations generate a deeper null towards the angles of interfering 
signals which is about (-61 dB for LA), (-60 dB for CA), and (-96 dB for PA). In the case of Rayleigh fading 
channel, the LA, CA, and PA configuration approximately converge after 20 iterations samples. The deeper 
null generated is (-62 dB for LA), (-48 dB for CA), and (-90 dB for PA) .The simulation results confirm that 
the performance of the RDES with PA has better main lobe with lower beam width but more side lobe inside 
of the appearance of great lobe effects, whereas CA is better in term of the number of side lobe than both LA 
and also PA antenna systems with fewer elements beside it's flexible angle range. 
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